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ABSTRACT 


The  plane  strain  solution  is  obtained  for  the  natural  vibrations  and  im¬ 
pulse  response  of  a  thin  circular  cylinder  containing  an  added  line  mass.  The 
solution  for  a  uniform  cylinder  is  derived  by  taking  the  added  mass  to  be  zero. 
Numerical  calculations  of  the  frequencies  and  mode  shapes  for  several  of  the 
lower  modes  are  presented  in  graphical  form  for  various  values  of  the  added 
mass.  The  general  impulse  response  solution. for  arbitrary  initial  conditions 
is  obtained  by  normal  mode  theory.  For  both  . the  natural  vibrations  and  im¬ 
pulse  response,  the  theory  is  found  to  bo  in  reasonnbloagroement  with 
available  experimental  results. 

In  a  particular  mode,  four  distinct  solution  statcsaro  found  to  exist: 
a  symmetrical  and  anti-symmetrical  branch  for  each  class  of  vibration,  flex¬ 
ural  and  extensional.  Noteworthy  features  revealed  by  this  investigation  are 
the  difference  in  frequency  and  mode  shape  of  each  solution  state  anil  the 
presence  of  couplingtbetwcen  the  flexural  and  extensional  classes,  particular¬ 
ly  noticeable  in  the  extensional  class  mode  shapes,  In  comparing  impulse  re¬ 
sponse  solutions  for  velocity  with  and  without  the  added  mass,  the  major  in-, 
fluenco  of  the  added  mass  is  found  to  be  an  increased  participation  of  the 
flexural  class  modes,  including  the  rigid  body  translation,  and  decreased 
participation  of  the  extensional  class  oscillatory  inodes. 


ADMINISTRATIVE  INFORMATION 

The  research  reported  herein  formed  a  thesis  submitted  to  and  approved  by  the  Graduate 
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degree  of  Doctor  of  Philosophy  in  Engineering  Mechanics, 

This  report  is  related  to  the  program  entitled  “Hull  Response  Inputs  to  Equipment,” 
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I.  INTRODUCTION  AND  LITERATURE. SURVEY 

The  vibrations  of  mass  loaded  cylindrical  shells  are  of  interest 
in  many  engineering  situations;  e.g.,  in  the  dynamic  behavior  of  hydro 
and  aerospace  vehicles  with  equipment  attached  to  their  cylindrical 
walls,  v/it'n  the  added  mass  considered  as  aq  imperfection,  such  a  study 
should  also  be  important  in  fundamental  studies  of  vibrations  of 
imperfect  bodies  of  revolution  inasmuch  as  the  effect  of  such  imperfec¬ 
tions  can  be  quite1  distinct  as  has  been  shown  by  several  experimental 
investigations  [1,2,3]*.  A  study  of  the  impulse  response  of  such 
structures  is  of  particular  interest  in  military  defense  applications, 
both  in  response  of  the-  shell  itself  and  of  attached  equipment  repre¬ 
sented  by  the  added  mass. 

The  first  formal  analytical  work  in  this  area  .apparently  was  by 
Den  riartog  [4]  in  1928.  In  studying  the  natural  vibrations  of  ring-like 
frames  of  electrical  machines  containing  lumped  masses,  he  derived 
approximate  fundamental  flexural  frequencies  for  these  frames  by  the 
Rayleigh  method.  In  1962,.  Palmer  [5]  investigated  the  variations  in  the 
natural  flexural  vibrations  without  extension-  due  to  a  lumped  mass  on  a 
ring  by  solving  the  classical  eigenvalue  problem.  A  significant  result 
of  that  study  was  the  difference  in  frequency  and  mode  shape  found 
between  the  symmetrical  and  anti-symmetrical  branches  of  each  mode. 
Michalopoulos  and  Muster  [6]  in  1966  obtained  the  natural  vibration 
solution  of  a  ring-stiffened,  mass-ioaded  cylindrical  shell  using  an 
energy  formulation  with  the  displacements  represented  approximately  by 

Numbers  in  brackets,  [  ],  refer  to  references  listed  on  page  64. 


finite  Fourier  series.  The  natural  vibrations  of  a  stiffened  pressur¬ 
ized  cylindrical  shell  with  an  attached  mass  were  considered  by  Ojalvo 
and  Newman  [7]  in  1967  by  an  approximate  energy  technique.  They 
determined  mode  shapes  for  the  shell  without  the  attached  mass  and  used 
these  as  assumed  mode  shapes  to  determine  the  natural  vibrations  for 
the  shell-mass  system.  Lee,  et  al,  [8]  in  1968  analyzed  the  natural 
vibrations  of  mass-loaded  ring  structures  by  two  approximate  methods, 
finite-element  and  transfer  matrix.  Using  essentially  the  same  approach 
as  Michalopoulos  and  Muster,  Lee,  et  al,  also  studied  the  natural 
vibrations  of  mass-loaded  shell  structures.  Although  the  natural 
vibrations  of  mass-loaded  beam  ,and  plate  structures  are  well  established, 
e.g. ,  Thornton  [9],  a  literature  search  has  failed  to  reveal  a  complete 
rigorous  solution  for  the  natural  vibrations  of  mass-loaded  ring  and 
cylindrical  shell  structures. 

Previous  solution  attempts  to  the  subject  problem  can,  with  one 
exception  [5],  be  characterized  by  the  approach  used  by  Ojalvo  and 
Newman  [7],  in  which  the  mode  shapes  for  the  uniform  shell  without 
added  mass  were  used  as  trial  functions  in  an  energy  formulation  for 
the  solution  of  the  complete  shell-mass  system.  A  rigorous  solution  as 
presented  herein  is  desirable  since  it  is  well  known  that  even  crude 
trial  functions  used  in  an  energy  formulation  can  give  quite  good  results 
for  the  natural  frequencies,  at  least  for  the  lower  modes,  but  mode 
shapes ,  and  response  parameters  computed  from  them  for  a  specified 
problem,  may  be  subject  to  question. 

The  problem  of  impulse  response  of  a  mass'.rloaded  cylindrical  shell 
has  been  treated  by  Michalopoulos  and  Muster  t'6].  They  used  modal 


analysis  to  formulate  the  response  of  the  structure  to  a  band  of  uniform 
infernal  pressure  applied  as  a  pulse  in  time.  However,  numerical  results 
were  not  presented.  With  the  above  exception,  the  impulse  response  of 
mass-loaded  shell  structures  apparently  has  not  been  treated  in  the 
literature.  The  simpler,  degenerate  problem  of  the  response  of  uniform 
thin  rings  and  cylinders  to  non-moving  impulse  loads' has,  however, 
received  some  attention.  Palmer  [10]  used  modal  analysis  to  obtain  the 
response  of  a  ring  to  a  side  pressure  pulse  represented  as  an  initial 
velocity  distribution.  In  that  study,,  the  effects  of  flexure  and  exten¬ 
sion  were  completely  uncoupled.  Humphreys  and  Winter  [11]  employed  the 
Laplace  transform  approach  to  solve  slightly  simplified  thin  shell 
equations  for  an  infinitely  long  cylinder  under  a  transverse  pressure 
pulse.  Sheng  [12]  solved  Donnell-type  shell  equations  for  the  response 
of  a  cylindrical  shell  to  a  semisinusoidal  surface  pulse.  Johnson  and 
Greif  [13]  used  two  different  numerical  methods  of  timewise  integration 
to  obtain  the  response  of  a  cylindrical  shell  to  a  blast  loading.  Their 
solution  was  completely  numerical  and  involved  expressing  Sanders' 
shell  equations  in  finite  difference  form. 
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II.  STATEMENT  'OF  THE  PROBLEM 

The  plane  strain  solution  for  the  natural  vibrations  and  impulse 
response  of  infinitely  long,  thiri\  circular  cylindrical  shells  containing 
an  added. concentrated'  line  mass'  will  be  obtained.  The  shell  material 
will  be  isotropic,  homogenous  and  the  analysis  will  employ  linearized, 
small  deformation,  thin  shell  theory,  neglecting  the  effects  of  shear 
deformation  and  rotatory  inertia,  as  formulated  by  Flugge  [14], 

The  midsurface  displacements  u,  v,  and  w  and  the  coordinate  system 
for  the- shell  are  defined  in  Figure  1.  The  line  mass,  M,  is  concen¬ 
trated  at  the  midsurface  and  is  parallel  to  the  z-axis.  Rotatory 
inertia  of  the  concentrated  mass  is  neglected.  Motion  will  occur  only 
in  planes  perpendicular  to  the  z-axis. 

The  solution  for  the  natural  frequencies,  and  mode  shapes  will  be 

obtained  in  the  classical  manner  by  describing  continuity  and  equilibrium 

conditions  at  the  mass.  The  influence  of  the  mass  will  be  ascertained 

by  comparison  with  the  vibrations  of  a  cylinder  without  added  mass. 

The  ex$£t-  solution  to  the  governing  equations  will  be  obtained  for  both 

extensional  and  flexural  classes  of  vibration  with  no  approximations 

made  on  either  class.  In  a  previous  solution -by  Palmer  [S']  only  the 

flexural  class  was  considered  in  an  approximate  manner,  i.e. ,  flexure 

3v 

with  the  midsurface  extension  everywhere  zero,  or  w  =  —  .  The 

30 

analytical  frequencies  will  be  compared  with  those  obtained  experimentally 
by  Lee  [8], 

The  solution  for  the  natural  vibrations  will  be  used  to  determine 
the  response  of  the  shell  to  a  transverse  impulse.  The  impulse  will  be 
specified  in  terms  of  initial  conditions  on  the  displacements  and 


velocities,  and  then  the  mode  constants  will  be  determined  through  the 
orthogonality  relation.  The  influence  of  the  mass  will  be  determined 
by  comparing  the  response  with  and  without  the  mass .  As  noted  in- the 
literature  review.  Palmer  [10-3  obtained  the  response  without  the  added 
mass  to  an  initial  velocity  distribution  by  modal  analysis.  Although 
both,  extensional  and-  flexural  classes  were  considered,  that  solution 
was  an  approximation  since  extension  was  assumed  to  occur  without  flexure 
and  flexure  without  extension.  Such  simplifications  permitted  a  solution 
to  be  obtained  in  a  relatively  straightforward-  manner.  The  impulse 
response  solution  will  be  compared  with  experimental  results  of  Navy 
tests,  .reported  in  [15]' and  [20],  and  Sandia  Corporation  tests  [16]', 


Figure  1  -  Coordinate  System 
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III.  NATURAL  VIBRATIONS 


The  general  solution  will  be  obtained  for  the  natural  frequencies 
and  normal  mode  shapes.  Numerical  calculations  of  the  frequencies  and 
mode  shapes  for  several'  of  the  lower  modes  will  be  presented  in  graphi¬ 
cal  form- and  compared  with  experimental  results  [8]. 


A.  Governing  equations 

The  complete  thin  cylindrical  shell  theory  equations  of  motion 
presented  by  Flugge  [14],  using  the  coordinate  system  of  Figure  1, 
adding  inertia  terms  and  omitting  load  terms  for  natural  vibrations,  are 


UZZ+ 


-v  .  \+v  K.  Fl-i>  ,  l-i>  "1  miia5  . 
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For  the  case  of  plane  strain,  u  =  0  and  (  )A~0  ,  these  three 
equations  reduce  to  two  in  the  planar  displacement  components  v  and  w  , 
Thus 


■7.. 

ma“v  , 

-dT+v«i> 


2- 

ma  w  , 

~dT  +  vo 


w0=0' 
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where 


and 
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Since  it  is  desired  to  treat  rings  as  well  as  stiffened  cylinders, 
■it  is  convenient  to  change  the  forv  .of  the  stiffness  parameters  Ds 
and  K  for  a  ring  of  unit  width  as!  follows : 


lVr?-EI,*EA 


and 


K_  =  J_  f  Eh3  1_  I 
L)sa2  l)sa2  1 2(  I  -  p“)  J  Aa2 


With  these  substitutions,  equations  (4)  and  (5)  can  be  written  in 
the  form 


hA 


-  IllilV  (Vqq  -  W<j)  =  0 


(6) 


lv\  FI 

mau  +  ~  (v0  -w)-~  (wmo  +  2w00  +  w)  =  0 
*1 


(7) 


It  is  noted  that  equations  (6)  and  (.7)  also  were  derived  by  Baron 
and  Bleich  [17]. 


B.  Solution  of  the  Governing  Equations 
Suitable  solutions  to  equations  (6)  and  (7)  for  harmonic  vibrations 
are ,  for  a  particular  mode , 

v  =  V(0)  eiwt  (8) 

and  w  =  W(0)  eit01  (9) 


where 


cos  iij.0  +  Bjj  sin  nj.0) 


(10) 


and 


"k)  (Aksin  nk(l-  Bkcosnk^ 


(11) 


Substituting  (8)  and  (9)  into  (6)  and  (7),  it  is  found  that 
equation  (6)  is  identically  satisfied,  whereas  equation  (7)  will  be 


satisfied  if  the  following  relation  is  in  turn  satisfied: 

C?-  |(n2-l)2Z  +  »2+  lJ.C  +  n2'(n2-  I)2  Z  =  0 

where 


!iA 


and 


Z=- 


Aa 


,2 


(12) 

(13) 

(14) 


The  n. ,  n„,  and  n  of  equations  (10)  and  (ll),  are  roots  of  equation 
Id  } 

2 

(12),  which  is  cubic  in  n  .  As  expected,  equation  (12),  the  frequency 
equation,  is  quadratic  in  the  frequency  parameter  C.  It  will  be  shown 
shortly  that  for  a  given  value  of  n  the  lesser  value  of  C  represents 
primarily  flexural  vibration  and  the  greater  value  represents  primarily 
extensional,  or  membrane,  vibration. 

C.  Consideration  of  the  Added  Mass 
As  the  solutions  (10)  and  (11)  contain  six  constants  of  integration, 
Ak  and  Bk  at  least  one  always  remaining  arbitrary  since  (6)  and  (7) 
are  homogeneous,  and  since  the  frequency  parameter  C  must  also  be 
determined,  then  six  relations  are  heeded  at  the  added  mass. 

From  continuity  conditions,  -the  following  equations  can  be  written 


M 


“  "|(J 

v  "n  =v(%i)-v(-0m)  =  o 

L 


w.l  M  =0 


H 

H 


-0. 


M  - 


=  0 


(15) ' 

(16) 
(17) 


rrom  equilibrium,  see  Figure  2,  summing  forces  in  the  radial 
direction  gives 

Mw(0m)  +  |N|  „  =0  (18) 


n°M 

-  N  =0 

L 
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and  in  the  tangential  direction  gives 


from  (15) 


from  (17) 


and  from  (18) 


where  R = 


M 


2rr:nn 


t 


B|.  nj*  sin  a^  =  0 


g  Bk  £tr|  (C  -  »£)  cos  ak  -<n£  sin  a  kj  =  0 


,  the  mass  ratio. 


(27) 


(28) 


(29) 

(30) 


and  «k  =  nk  0^.  k  =  1 , 2,  or  3. 


(31) 


Without  any  loss,  of  generality,  0^  will  be  taken  as  jt  for  the 
remainder  of  this  dissertation. 

D.  Natural  Frequencies 

Since  (24)  through  (29)  are  homogeneous,  the  determinant  of  the 
coefficients  of  the  Ak  and  Bk  terms  must  vanish  in  order  for  non¬ 
trivial  values  of  these  constants  to  exist. 

The  sixth  order  determinant  formed  from  (24)  through  (29)  is 


l)  = 


0 


(32) 


12 


I  V' 


where 


•>A" 


(C  ny)  sin  «|  (C  n;).®^  (C  up  sin  <13 


•>  ^ 
n Y  (C  -  up  sin  a  j 


i  _  ■) 

no  (C  im)  snt«i 


■> 

n-pf  up  sin  (13 


n  |  (nUCcosaj  +ii|  sin  up  n->  (nRC  cosa-,  4  n->Miia-,)  113  (nRC  cosaj  +  113  sin  «p 


(33) 


and 


%= 


n  |  sin  a  | 


n  |  sin  a  | 


-.111  Sill  tti 


n->  sin  a-> 


11^  sin  a  t 


«UrpC  -  up  cos  a | 


ttR^-CG*  -  11^)  cos  a-> 
-  Ht  sin  a-. 


113  sin  <13 


n‘3  sin  «3 


-i>C 


K~(C  -  up  cos  <13 
-  113  sin  03 


(34) 


Equation  (32)  is  expressed  as  a  product  of  two  third-order  determi¬ 
nants,  one  containing  only  coefficients  of  Ak  and  the  other  containing 
only  coefficients  of  B,.  ,  Thus 


or 


D  =  DaDb  =  0  (35) 

Equation,  ( 35)  may  be  satisfied  in  three  ways 


da-=0 

(36) 

db  =  o 

(37) 

Da  =  D[3=0 

(38) 

In  the  first  of  these,  (.36).,  nontrivial  values  are  possible  for  the 
Aj.  constants  since  DA  =  0  .  But  if  DA  =  0  and  D^O  ,  then 

all  the  Bk  constants  will  be  zero,,  since  by  Cramer’s  rule  nontrivial 


values  cannot  exist.  The  reverse  is  true  for  the  second  of  these,  (37). 
Referring  to  the  expression  for  the  radial  displacement  function,  (11), 
we  find  that  the  Aj.  constants  are  associated  with  sine  functions  which 
are  antisymmetrical ,  whereas  the  constants  are  associated  with 

symmetrical  cosine  functions.  Condition  (36)  then  would  result  in  a 
solution  containing  only  antisymmetrical  terms  and  the  solution  from  (37) 
would  contain  only  symmetrical  terms. 

in  either  case,  (36)  or  (37),  both  the  governing  equations  and. 
the  continuity  and  equilibrium  conditions  at  the  mass  will  be  satisfied. 

The  frequency  equations  may  be  found  by  evaluating  DA  and  D[}  • 

from  ( 33) 

O  O  } 

D  ^  =  n yn->  (C  ~  ny)(C-  115)  sin  a  |  sin  (nRC  cos  a->  +  n^  sin  a-,) 

0  0  0 

+  11:5113  (C  -  ny)  (C  -  115)  sin  «|  sin  «->  (mRC  COSCI3  +  113  sin  a-j) 

+  n?n.(C-  115)  (C  -  ni).sina->  sin  a^(^RC  cosn.  +  11,  sin  a,) 

■  ;  ;  “■  (39) 

-  n |  ns>  (C  -  iii)  (C  -  np  sin  a-,  sin  a  3  (mRC  cos  (tj  +  n  |  sin  aj ) 

O  0  -> 

-  n^n-f  (C  -  ny)  (C  -  np  sin  <i|  sin  (13  (;?RC  cos  a-,  +  n7  sin  a-,) 

0  o  o 

-  n<jny(C-  np(C-  n^)  sin  a  j,  sin  a->  (ttRC  cos  U3  +  113  sin  <13) 

Setting  =  0  the  frequency  equation  for  antisymmetrical  vibra¬ 
tions  is  obtained.  Solving  this  expression  for  the  mass  ratio  and 
letting  R  =  Ra  for  DA  =  0  ?^ves 


Ra  =  »? 


A  ~  «f  (C  -  up  (C  -  !i5)(n5  -  n?)  +  115  (C  -  nf)  (C  - 115)  (113  -  ny) 


+  115  (C  -  ny)  (C  -  115)  (ny  -  npj*  nC  jn |  (C  -  115)  (C  -  n j)  (115  -  115) 

+  «»2  (C  -  ny)  (C  -  115)  (ny  -  np  -r^-^  +  113  (C  -  ny)  (C  -  np  (115  -  ny)  yr— ] 


(40) 


//. 


By  the  same  procedure,  we  find  from  (37) 

dB  =  nin2n3  [n2n3  <n2  '  n3^ +  n?n2  (n?  "  n2> +  n?n3  (n3  “  n  1^] 

T  „  cosa3  to  2  cosal 

+  jtr|  jnji^ (n2~  "D  (C  '  n3>  IhT^  +  n2n3  (n3  "  n2>  <C  "  n I ^ 

* .3(»?-"|)(C-i)CS]  <41) 

from  which  the  mass  ratio  for  symmetrical  vibrations,  Rg  ,  is 
rb  =  n1n2n3  (115  -  115)  +  »,j "2  ("2  “  nl> +  “Tn3  ("?  '  n5)] 

-  4  [n2n3  ("3  '  "2>  (C  "  n  I  ^  +  "  1  n3  (n  I  -  "3>rC-  "2*^ 

-+;n1n2(4-?)(C-,i)^.  ^ 

In  general,  both  DA  and,  DB  will  not  vanish  for  the  same  nk  ,  C 
and  R  values  since  expressions  (40)  and  (42)  are  not  the  same.  New 
referring  to  equations  (12).,  (40),  and.  (42)  ,  it  is  apparent  that  for 
natural  vibrations  with  an  added  mass  four  solutions  are  possible;  one 
consisting  of  symmetrical  vibrations  and  another  of  antisymmetrical 
vibrations  for  each  of  the  two  classes,  flexural  and  extensional,  An 
illustration  of  this  is  provided  in  the  next  section  where,  numerical 
solutions  are  obtained. 

The  frequency  equation  for  the  case  of  a  uniform  cylinder  or  ring 
may  be  obtained  by  letting  the  mass  ratio  become  zero  in  (40)  and  (42). 
From  (40),  multiplying  through  by  sin  a,  sin  a2  sin  a3  t 
jn|  (C  -  nJj)  (C  -  ny)  (115  -  115)  +  w\  (C  -  nj)  (C  -  n|)  (n|  -  ny) 

+  n|(C  -  ny)  (C  -n^)  (ny-n^jj  sin  a,  sin  a,  sin  a3  =  0  (4- 

and  from  (42),  again  multiplying  by  sin  q,  -sin  a2  sin  a3  * 

n2n3  (n3  ~  11 2^  +  nIn2  (n2  "  nP 

+  n^n3  (nj  -  n^jij  sin  a|  sin  a2  sin  a3  =  0 


(44) 


,  if  one  of  the  is  an  integer, 


In  this  case,  with  0^  =  n 
then  (43)  or  (44)  will  be  satisfied.  Actually,  for  this  situation,  both 
(43)  and  (44)  are  identically  zero,  which  is  condition  (38).  Then 
equation  (12)  suffices  as  the  frequency  equation,  with  one  of  the  n^  an 

integer. 

L.  Normal  Mode  Shapes 

The  antisymmetrical  branch  mode  constants  Aj.  may  be  found  by 
solving  (24)  and  (25)  simultaneously;  the  solution  is 


and 


0  0  o 

At  (nj  -  ny)  (C  -  n]")  sin  a( 

"  ~  0  0  0 

A  |  (115  -  iip  (C  -  n  t)  sin  «t 


00  o 
A^  (nf  -  »t)  (G  -  np  sin 

A  |  (115  -  up  (C  -  nj)  sin  <13 


(45) 


(46) 


and  the  symmetrical  constants  Bk  from  (27)  and  (28)  are 

0  o 

B-,  11 1  (115-  up  sin  a | 

B|  iit  (n?  -  113)  sin  a0 


(47) 


and 


B3  11 1  (nj-  -  115)  sin  <i| 

—  -  ^  ^ 

Bj  113  (»t  -  up  sin  a3 


(48) 


where  and  3^  are  arbitrary.  From  (10)  And  (11)  the  tangential  and 

radial  displacement  functions  for  the  antisymmetrical  branch  axxj 

f  At  a3  1 

n  |  cos 11  \0  +  ^  ib  cos  Dt0  +  —  113  cos  113OJ  (49), 


and 


V  =  A , 


W  =  A 


T  At  T  ^3  - 

(C  -  np  sin  ii|0  +^(C  up  sin  +  -^ — (C  - -up  sin  ivjflj 


(50) 


Ifi. 
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Jh\  l"s',n  ("  I  is  ‘  n3is) n  sin(,,lis  +  "3is^ 

Ylll“\B1jnlisn3is[  'Mis  ’  n3is  '  "  I  is +  n3is 


siu  <n2is "  n3is>^  <fn  (n2is  +  n3is>ff 
n2is  "  n3is  n2is  +  n3is 


F,  Numerical  Solutions 


To  illustrate  the  influence  of  the  added  mass,  frequencies  and  mode 
shapes  were  calculated  for  several  of  the  lower  modes  by  the  following 
procedure : 

1.  Influence  of  Mass  on  frequencies 

a.  After  selecting  a  value  of  the  stiffness  parameter  Z, 


values  of  n^  wei’e  chosen  and  (12)  was  solved  for  the  flexural  frequency 


parameters  and  the  extensional  frequency  parameters  Cg.  Then, for 
each  frequency  parameter,  and  Cg,  (12)  was  solved  again  for  the  modal 
parameters  n^  and  n^  for  the  flexural  class,  and  n^  and  n3g  for  the 


extensional  class.  Integer  values  of  n^  will  give  solutions  for  the 


problem  without  the  added  mass. 


b.  The  n^  and  C  values  were  then  inserted  into  equations  (40) 


and  (42)  to  determine  values  of.  the  mass  ratio,  R.  Note  that  for  each 

n.  chosen,  four  values  of  R  will  be  obtained,  R_  ,  R.  ,  R  and  R  , 

1  icl  XS  €3/  GS 

where  f  corresponds  to  flexure,  e  to  extension,  a  to  antisymmetric  and 
s  to  symmetric. 

c.  Frequency  ratios  of  the  structure  with  added  mass  to  that 
without  added  mass  were  obtained  from  equation  (13): 

.1/2 


_  NV 

wm  \^m  / 


(89) 


For  several  individual  lower  modes ,  the  above  parameters  were 
calculated  and  are  given  in  Tables  1  and  2.  The  relationship  between 
frequency  and  mass  is  also  shown  graphically  in  Figures  3  through  6. 
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TABLE  1  -  NATURAL  FREQUENCIES,  EXTENSIONAL  CLASS,  7,  =  0,0001 


n, 

C 

R 

a  a 

R 

AC 

ZEROTlt 

MODE 

0.0 

1.000100 

9.9247771 

10.074785 

1.0 

0.0 

0.0 

0.093543i 

0.991352 

9.9026671 

10.053005 

0.996 

— 

0.10 

0.096365i 

0.990817 

9 .9013081 

10.051667 

0.995 

— 

1.0 

0.096652i 

0.990761 

9.9011681 

10.051528 

0.995 

— 

10.0 

TIRST 

MODE 

1.0 

2.0 

11.8288731 

11.955009 

1.0 

0.0 

0.0 

0.920444 

1.847213 

11.5936511 

11.722319 

0.962 

0.086 

0.10 

0.909462 

1.827122 

11. 561633i 

11.690650 

0.956 

0.10 

0.178 

0.895314 

1.801590 

11. 520569i 

11.650043 

0.950 

0.118 

1.0 

0.892210 

1.796041 

11. 5il586i 

11.461161 

0.948 

0.123 

10.0 

0.685508 

1.469940 

10.9425731 

11.078805 

0.858 

1.0 

—  — 

0.668459 

1.446858 

10.898814i 

11.035586 

0.852 

10.0 

— 

SECOND 

MODE 

2.0 

5.000180 

14. 90 362 3i 

15.003939 

1.0 

0.0 

0.0 

1.923135 

4.698604 

14, 672060 i 

14.773947 

0.970 

0.04 

0.10 

1.905064 

4.629420 

14.617369i 

14.719634 

0.962 

0.051 

1.0 

1.902883 

4.621112 

14.6107601 

14.713071 

0.961 

0.052 

10.0 

1.832289 

4.357412 

14.3961911 

14.500016 

0.934 

0.10 

-- 

1.673354 

3.800200 

13,908495i 

14.015931 

0.872 

.1.0 

— 

1.671937 

3.795457 

13.9041191 

14.011589 

0.871 

10.0 

— 
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TABLE  2  -  NATURAL  FREQUENCIES,  I LEXURAL  CLASS,  Z  =  0.0001 


nl 

C 

n2 

n3 

“Nl 

“m 

R, 

fa 

R, 

fs 

2.0 

0.000720 

SECOND 
(n2  =  a+bi, 

a 

0.4132*1 

MODE 

n3  =  a-bi) 

b 

1.081867 

1.0 

0.0 

0.0 

1.988584 

0.000697 

0.418348 

1.073266 

0.984 

0.10 

0.021 

1.956736 

0.000634 

0.431774 

1.049056 

0.939 

1.0 

0.089 

1.952292 

0.000626 

0.433560 

1.045649 

0.933 

1.43 

0.10 

1.939996 

0.000603 

0.438395 

1.036185 

0.915 

10.0 

0.132 

1.798911 

0.000382 

0.483953 

0.923071 

0.728 

— 

1.0 

1.705399 

0.000271 

0.505549 

0.842400 

0.614 

10.0 

3.0 

0.005760 

THIRD 

1.0368411 

MODE 

2 ,432945i 

1.0 

0.0 

0.0 

2.989798 

0.005668 

1.0376341 

2.420029i 

0.992' 

0.10 

0.009 

2.963754 

0.005439 

1.039746i 

2.386907i 

0.972 

1.0 

0.035 

2.951385 

0.005333 

1.040797i 

2.371093i 

0.962 

10.0 

0.048 

2.911696 

0.005002 

1.044397i 

2.319959i 

0.932 

— 

0.10 

2.737457 

0.003720 

1.0658711 

2 .0865931 

0.804 

-r 

1.0 

2.680194 

0.003356 

1.0759711 

2.0055851 

0.763 

— 

10.0 

Variation  of  Frequency  with  Mass,,  Zeroth  and  First  Modes 
Extensional  Class,  Z  =  0.0001 


SHE 


Figure  5  -  Variation  of  Frequency  with  Mass,  Second  Mode,  Flexural 

Class,  Z  =  0.0001 


Class,  Z  =  0.0001 


As  noted  on  the  figures,  the  stiffness  parameter  7.  was  taken  to  be 

lO-'*.  This  particular  value  of  7.  represents  a  ring  of  unit  width,  12 

inches  in  diameter,  with  a  thickness  of  0.21  inches.  Calculations  using 

—2  -6 

*-th«r  v ilues  of  7.  ranging  from  10  to  10  showed  less  than  10%  differ- 
*’nro  with  those  represented  by  the  figures,  with  the  largest  differences 
occurring  in  the  extensional  class  modes  and  virtually  no  difference  in 
"he  ‘ iexural  class  particularly  for  smaller  Z  values. 

•  n*. Iuo"*e  of  Mass  on  Mode  Shapes 

equations  (49)  through  (52)  were  used  to  determine  the  normal¬ 
ised  mode  shapes  for  desired  values  of  the  mass  ratio,  R.  The  parame¬ 
ters,  n^_  and  C,  for  particular  values  of  R  were  found  by  iteration  on 
the  process  described  previously  in  l.a.  and  l.b. 

The  normalized  displacement  function  »  is  shown  in  Figures  7 
and  3  for  0  between  0  and  180-  degrees. 

G.  Comparison  with  Experiment 

Lee,  et  al,  [8]  conducted  an  experimental'  study,  -to  determine  the 
effects  of  added  mass  on  the  vibrations  of:  ring  structures.  By  reson¬ 
ance  testing,  they  found  flexural  symmetric  natural  frequencies  and 
mode  shapes  for  several  values  of  the  mass  ratio,  R.  Their  frequency 
results  are  shown  on  Figures  5  and  6  for  comparison  with  present  theory. 
The  experimental  results  show  the  same  general  trend  as  theory  for 
increasing  P,  but  the  experimental  frequency  ratio  is,  on  the  average, 
about  10 %  higher  than  theory.  For  a  variety  of  reasons,  precise  agree¬ 
ment  between  theory  and  experiment  generally  cannot  be  expected, 
especially  for  vibratory  and  transient  phenomena.  In  this  case  the 
-  rincioal  cause  of  disagreement  can  be  attributed  to  the  fact  thv.t  in 

os 


SECOND  MODE,  SYMMETRICAL  BRANCH 


Figure  7  -  Normalized  Mode  Shapes,  Extensional  Class 

■  nnm 


the  conduct  of  the  experiments  the  ring  was  driven  eccentrically  with 
respect  to  the  attached  mass,  thus  permitting  some  participation  of  the 
antisymmetric  branches.  As  shown  by  present  theory,  this  would  account 
for  somewhat  higher  frequency  ratios  than  those  that  would  result  from 
pure  symmetric  excitation.  Experimental  mode  shapes  are  shown  in 
Figures  9  and  10  and  show  the  same  general  trend  of  reduced  amplitude 
at  the  mass  for  increasing  R  as  the  theory. 
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R  =  1.58 


R  =  3.52 


Figure  9  -  Experimental  ['8j  Mode  Shapes,  Second  Mode 
Flexural  Class,  Symmetrical  Branch,  Z  =  2.6  x  10_s 


Figure  10  -  Experimental  [8]  Mode  Shapes,  Third  Mode, 
Flexural  Class,  Symmetrical  Branch,  Z  =  2.6  x  10"6 


11.  Discussion 


The  influence  of  the  added  mass  on  the  natural  vibrations  will  be 
discussed  in  terms  of  the  solution  without  added  mass. 

1.  Influence  of  Mass  on  frequencies 

An  obvious  feature  here  is  the  difference  in  frequency  between 
the  symmetrical  and  antisymroetrical  brancnes  of  a  normal  mode,  as  shown 
in  figures  3  through  6.  This  holds  for  both  classes  of  vibration, 
extensional  and  flexural,  for  no  added  mass ,  R  =  0,  both  branches  of  a 
particular  class  have  the  same  frequency,  but  as  the  added  mass  increas¬ 
es,  the  difference  in  frequency  also  increases.  Note  that  the  zeroth, 
or  breathing,  mode  has  only  the  symmetrical  branch.  This  is  an  expected 
result  since  for  the  case  of  zero  added  mass,  only  the  symmetrical 
branch  exists,  as  shown  by  equations  (53'>  through  (56). 

for  the  extensional  class  modes,  as  the  mass  ratio  increases, 
the  reduction  in  frequency  is  greater  in  the  antisymmetrical  branch 
than  in  the  symmetrical  branch.  The  reverse  is  true  for  the  flexural 
class. 

With  the  exception  of  the  zeroth  mode,  it  is  noted  that  for 
either  class,  the  reduction  in  frequency  is  less  for  both  branches  for 
the  higher  modes  than  for  the  lower.  The  zeroth  mode  shows  very  little 
reduction  in  frequency  regardless  of  the  amount  of  added  mass. 

2 .  Influence  of  Mass  on  “ode  Chapes 

Figures  7  and  8  indicate  that  the  displacement  of  the  added 
mass  decreases  as  the  mass  ratio  increases.  A  large  added  mass  tends 
to  become  a  node.  In  the  flexural  class,  the  mass  has  appreciably  less 
influence  on  the  antisymmetrical  shapes  than  on  the  symmetrical,  as 

■u 


shown  by  Figure  8.  The  effect  of  added  mass  on  the  extensional  class 
mode  shapes  is  striking.  Figure  7  shows  the  presence  of  coupling  with 
higher  flexural  class  modes.  F.ven  for  a  very  small  value  of  the  mass 
ratio  (R  =  0.01)  this  coupling  is  quite  severe,  particularly  for  the 
symmetrical  branches.  Observe  also  that  even  with  this  coupling  the 
true  character  of  the  extensional  shapes  is  maintained  since  the  flexure 
is  superimposed  on  the  primary  extension.  This  behavior  is  especially 
noticeable  in  the  zeroth  mode  shape  shown  in  Figure  7. 

For  the  uniform  cylinder  without  added  mass,  it  is  interesting  to 
note  from  equations  (54)  and  (56)  that  the  only  distinction  in  shape 
between  the  symmetrical  and  anti-symmetrical  branches  of  a  particular 
mode  is  that  the  axis  of  symmetry  of  the  anti-symmetrical  branch  shape 
is  rotated  ■—  radians  from  the  axis  of  symmetry  of  the  symmetrical 
branch  shape.  As  shown  in  Tables  1  and  2,  for  R  =  0,  both  branches 
of  the  same  class  of  vibration  have  the  same  natural  frequency.  Such 
degeneracy  does  not  exist  when  the  added  mass  is  present.  Not  only  are 
the  anti-symmetrical  branch  shapes  completely  without  symmetry  and 
different  from  those  of  the  symmetrical  branch,  but  the  frequencies 
associated  with  each  branch  are  different  as  well. 

The  zeroth  and  first  flexural  class  modes  are  non-oscillatory , 
rigid  body  modes.  They  will  be  treated  in  the  next  Chapter, 


IV.  IMPULSE  RL3I0H3K 


The  general  impulse  response  solution  will  be  obtained  by  modal 
analysis  for  arbitrary  spatial  variation  of  the  initial  conditions. 
This  general  solution  will  then  be  specialized  to  a  particular  form  of 
the  initial  conditions  for  comparison  with  experimental  results.  The 
impuLse  response  of  a  uniform  cylinder  (R  =  0)  will  be  obtained  as  a 
degenerate  case  of  the  general  impulse  response  solution, 

A.  General  Solution 

As  was  shown  in  Chapter  Ill,  for  each  mode  of  vibration  (zeroth, 

first,  etc.),  four  solutions  are  possible,  symmetric  and  antisymmetric 

branches  for  each  of  the  two  classes,  flexural  and  extensional,  Since 

each  of  the  four  is  a  solution  and  since  the  problem  is  linear,  then  a 

linear  combination  of  the  four  solutions  is  also  a  solution.  There- 

till 

lore,  the  displacement  components  for  the  l  mode  can  be  written  as 


vi  =veai  +  vesi  +  vlai  +  vls. 


w.Bwcai  +  wc»  +  wfai  +  wfd  (91) 

It  is  convenient  to  express  the  spatial  and  temporal  variables  of 

each  of  the  components  of  ei^uations  (90)  and  (91)  as  follows 


vea> 

"  Veai 

(c,  sin  <Jcait  f  (lj  cos  o;cait) 

(92) 

vesi 

=  Ves. 

(Cj  sin  oj0sjl  +  dj  cos  wcsjt) 

(93) 

'  f'.n 

=  V»;u 

(a.  sin  Cu.-yjt  +  b,  cos  oj|ajt) 

(94) 

V. 

"  Vl\. 

(.ij  sin  W|Sjl  +  bj  cos  Wfsjt) 

(95) 

weai  =  Weai  <Ci sin  weai<  +  di  cos  wcail>  (%) 

wcsi  "Wesi(°i siu  ^osi*  +  t!i  cos  tocsjt)  (97) 

wfai  =  Wfai  sin  wf;li(  +  b;  coscofajt)  (98) 

wfsi  =  Wfsi(:li  sin  wlsit  +  bj  coswrsjt)  (99) 


where  V  .  ,  W  .  ,  etc.,  are  the  mode  shapes  obtained  previously,  and 
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a^,  a^,  b etc.,  are  constants  to  be  obtained  from  the  initial 
conditions.  Summing  all  the  modes  gives  the  complete  solution  which  is 


and 


v  =£\;  (100) 

1=0 

w*£Wi  (101) 

i=0 


For  given  initial  conditions  on  the  displacement  and  velocity,  the 
constants  a^,  a^,  b^,  etc.,  can  be  obtained  through  the  orthogonality 
relation.  Let  the  initial  conditions  be  stated  as 


v(fl,0)  =v 
w(0 ,0)  =  w 

(102) 

v(0.O)  =v 
w(0,O)  =  w 


Then,  from  equations  (100)  and  (101),  with  (102)  we  have 
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(103) 


.38 


(Ill) 


^  ’/i(7T,0)  =v(7l) 
i^O 


^"^W^TT.O)  =  W(7l) 
i=0 


(112) 


y~^Vj(?r.O)  =v(w) 

i=0 


(113) 


^^Wj(jT,0)  =  w(fl) 
i=0 


(114) 


Multiplying  equations  (111)  and  (113)  by  MVfsj  and  (112)  and 


.(nr) 


(114)  by  M\V|y!j  gives 


^Mv.(n,0)Vg  =  Mv(fl)V§ 
i»0 

Yw*  Vfsj  =  Mv(7r)  Vfsj 


i=0 


^Mw^,0)wg  =  Mw(,)\V(g 

i=0 

oo 

^  ^NlWj(ff.O)  Wfjj  =  Mw(ff)  Wfsj 


i=0 


(115) 


(116) 


(117) 


(118) 


Adding  equations  (107),  (109),  (115)  and  (117)  results  in 
n  [Vj(0 .0)  Vfsj  +  Wj(0,O)  Wfsj]  madO  +  M  [vj(0)  Vfsj  .  +  Wj(0)  Wfj  ^  | 


i=0  W-v 


=  jf  [vV ,sj  +  wW,sj]  madO  +  M  [vVfsj  +  w\Vf J  ^ 


(119) 


Using  equations  (90)  through  (’J9),  and  the  orthogonality  condition, 


vV  •  +  w\V  ■ 
v  vesi  vv esi 


r  [*V6it*Wfal]  "tR[™firit*wrai] 


■£["***•*]»**[ 


Wfai+^Vfai 


U  [*-♦"-]  "tR[* 


vV.  ■  +  wW  • 
eai  wneai 


vV  ■  +  w\V 
cai  o 


Equations  (123),  (124),  (127),  and  (128)  associated  with  the 
extensional  class  are  valid  for  any  mode  number,  i,  equal  to  or  greater 
than  zero.  Equations  (121),  (122),  (125)  and  (126)  associated  with  the 
flexural  class  are  valid  only  for  mode  numbers  equal  to  or  greater 
than  two.  As  shown  previously  in  Chapter  III,  harmonic  solutions  are 
not  possible  for  the  zeroth  and  first  modes  of  the  flexural  class. 
Therefore  these  must  be  treated  separately.  Previous  analytical  studies 
[10,19]  of  the  flexural  response  of  uniform  rings  have  shown  that  the 
zeroth  and  first  flexural  clas*i  modes  are  associated  with  rigid  body 
motions.  The  solution  for  the  zeroth  mode  can  be  taken  as  follows 


vfO  =  fO  +  hOt  and  wf0  =  80  +  k0t 
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Substituting  U?yj  into  the  governing  equations  (6)  and  (7),  it  is 
found  that  g  =  -  0  and  f  an<l  h  ^  dre  arbitrary,  which  represents  a 

rigid  body  rotation.  For  the  first  mode,  the  solution  is  taken  as 


\|l  U(  t  It | n  \in/J  and  \\j|  =  (g|  a  kjl)  cosO 


(130) 


upon  substituting  (130)  into  equations  (6)  and  (7),  it  is  found 
that  the  governing  equations  are  satisfied  by  equations  (130)  if 
fl  '  ‘"l  anc^  =  ^l*  '*ote  that  this  solution  represents  a  rigid  body 
translation.  Tne  non-harmonic  or  zero  frequency  solutions  for  the 
zeroth  and  first  flexural  class  modes  given  by  equations  (129)  and  (130) 
also  satisfy  the  equilibrium  and  continuity  conditions,  equations  (15) 
through  (20),  thus  are  valid  solutions  to  the  given  problem.  The 
orthogonality  relation,  equation  (57),  also  is  applicable  for  these 
solutions . 

Froceeding  in  the  sane  manner  as  before,  the  constants  f  ,  h  ,  f, 

o’  o*  1 

and  h^  are  found  to  be 


s/: 


vdtf  +  Rv(jt) 


I  a  R 


4  jr* 


dfl  +  Rv(tt) 


h» r  I  +  R 


ifj 

if  ■ 


v  sin  0  t  w  cosfl]  <i 0  -  Rw(m) 


I  t  R 


v  sin  0  +  w  cos  0]  d0  -  Rw(m) 


I  +R 
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(131) 


(132) 


(133) 


(134) 


Thus,  with  all  the  constants  determined  from  equations  (121) 
through  (128)  and  (131)  and  (134),  the  complete  general  solution  is 


«•  m 

v  =  fQ  +  hQt  +  (f,  -Hi ,  t)  sin  0  £Vfi  £)< 

i-2  i«0 


w  =  (f|  +  h 1 1)  cos  0  + 


'X/Vfi  +Swci 


>2  F*0 


where 


vfi  =  vfai +  vfsi 


v  *  =  V  .  +  V  • 

vci  cai  csi 


Wfi  =  wrai+wfsi 


W  -  =  W  .  4'  W  * 

ei  cai  csi 


The  impulse  response  solution  for  a  uniform  cylinder  without  added 
mass  has  the  same  form  as  above  and  may  be  derived  from  the  above  gen¬ 
eral  solution  by  setting  R  =  0  and  using  the  mode  shapes  given  by 
equations  (53)  through  (56). 

B.  Solution  for  Specified  Initial  Conditions 
For  later  comparison  with  experimental  results,  the  impulse 
response  solution  for  a  particular  set  of  initial  conditions  will  now 
be  obtained.  Let  the  initial  conditions  be  of  the  following  form: 


v  = v  =w  =  0 


for  all  0 


w  =  -WqCosO  ,  for  7T  >0>-~ 

and  for~>0 ■>-  7r 
w  =  0  ,  for— 
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where 


frQ-W  """  ,  nktsi> "  Ml'-WS 

211-W  *  2(i*'W  ‘  2('-W  ‘ 

sin(l -nkesi)7t  sin(l +  n1.csj)n  sin  (1  -  nkcsj)§  sin  (I  +  nkcsi)f 

ki  -  2d-nkesir  +  2d  +  n^csj)  '  2(1  -n^F  '  2(fF^F 
k  =  1.2,3 


OQ  CO 

»'=lMtcos0  +  £  ajW^i  sin  wfsjt  +  £  CjWesi  sin  wesit  ( 142) 

i=2  i =0 

In  the  above  treatment,  it  is  emphasized  that  the  mode  shapes  and  Vd 
are  normalized  in  accordance  with  equation  (59), 

For  the  case  of  a  uniform  cylinder,  R  =  0,  the  form  of  solution 
given  by  equations  (141)  and  (142)  remains  the  same;  however,  the 
expressions  for  the  constants  h^,  a^  and  c^  are  simplified.  Equation 
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(  L3u)  becomes 


It 


I 


(143) 


equation  (139)  becomes 


iTcO 


fell 


I  sir  In 


n  Z* 


(144) 


and  equation  (140)  becomes 


wOnlcsn 

"^esn 


(,r-Cosn)Hl.i-  n  55  0 


(145) 


Note  that  the  mode  number  now  is  the  same  as  the  integer  n. 

Other  response  parameters,  such  as  velocity  and  stress,  can  be 
obtained  by  appropriate  temporal  and  spatial  differentiation  of  the 
displacements  given  by  equations  (141)  and  (142). 

C.  Comparison  with  Experimental  Results 
Several  specific  numerical  solutions  will  now  be  obtained  for 
comparison  with  available  experimental  results. 

1.  Cylinder  with  Added  Mass 

The  cylinder  is  impulsively  loaded  by  a  shockwave  from  a  near¬ 
by  underwater  explosion.  The  pertinent  parameters  for  this  problem  are 
u  =  20  x  10 k  (psi) 
v  =  0.3 
a  =  8.65  ( f t ) 

in  =  0.0502  (lb  ?-Q-c~  ■) 
in 

Z  =  0.00194 


It  is  desired  to  compare  with  Navy  experimental  radial  velocity 
measurements  recorded  at  the  location  of  the  added  mass,  o  -  n 
Differentiating  equation  (142)  with  respect  to  time,  we  obtain  the 
following  expression  for  the  non-dimensional  radial  velocity  at  0-v 


w(jl,  t) 

*0 


+  R 


I  +  R 


Ea;-w 


(H6) 


where 


w<*»  CO,- 

"I  si  lM 


Wn 


and 


"esi  wcm 


w0 

The  amplitude  constants  and  c|  are  listed  in  Table  3  through 
eight  modes  (i  =  8).  A  measure  of  convergence  is  indicated  by  the 
values  shown  in  Table  3  since  both  a^  and  c^  become  very  small  for  the 
higher  modes.  In  fact,  these  amplitude  constants  are  insignificant  for 
i  >  4.  Figure  11  shows  the  solution,  equation  (146),  summed  through 
eight  modes.  Shown  also  in  Figure  11  is  the  experimental  velocity 
history  from  reference  [20 ].  The  high  frequency  oscillations  on  the 
experimental  histories  are  due  primarily  to  the  response  of  the  velocity 
meter  elements.  The  fundamental  natural  frequencies  of  the  elements 
were  measured  and  found  to  correspond  with  the  superimposed  high- 
frequency  oscillations  of  the  experimental  records.  This  is  discussed 
further  in  reference  [21],  together  with  other  details  of  the  Kavy 
experimental  program. 


TABLE  3  -  VELOCITY  MODE  AMPLITUDES ,  0  =  n  ,  Z  *  0.00194 


0 

0 

0.3183 

0 

0.2112 

1 

0.25 

0.25 

0.3644 

0.0917 

2 

0.1696 

0.0424 

0.2534 

0.0067 

3 

0 

0 

0.0829 

-0.0004 

4 

-0.0398 

-0.0026 

0.0652 

-0.0013 

5 

0 

0, 

0 

-0;0008 

6 

0.0176 

0.0006 

-0,0289 

-0.0005 

7 

■o 

0 

-0.0064 

0 

8 

-0.0099 

-0.0002 

-6.0084 

0 

TIME(MSEC) 

Figure  11  -  Theoretical  and  Experimental  Velocity  Histories 

6  =  ti.  Z  =  0.00194 


2.  Cylinder  without  Added. Hass 

The  parameters  of  the  previous  problem  apply  here,  except  that 


R  =  0.  The  solution  in  *his  case  is 


w(rr.l)  t  ,  \  '  i  ,  ,  > 

T-  =  4  +  Zj  11,1 005  WfS"  2L'C» 

W0  II  =  2  11  =  0 


C0SC*W 


where 


2(-l)“  (n“  -Ct-sn)“ 

- - - - - - - f  n  -  2,  4,  6,  -  -  - 

n(l  -ir)  [(n“..-Cfsn)-  +  ir] 


an=  ,,  :  o 


aj,  =  0 ,  n  =  3.5,7.--. 


,  1 
cl=4 


...  2(-l)V.-Cctn)2 

S  itCI -n2)  l(n2-Ccsn)2  +  n2|  ’  " 


c'  =0.  n  =  3,  5,  7,  -  -  • 


(147) 


The  above  amplitude  constants  are  listed  in  Table.  3  for  eight  modes 

(n  =  8).  Both  of  the  above  expressions  for  a'  and  o'  rapidly  become 

n 

very  small  for  large  n  as  is  shown  by  the.  values  in  Table  3.  The 
velocity  history  computed  from-  equation  (147)  summed  through  eight  modes 
is  shown  in'  Figure  11,  together  with  the  experimental  velocity  history 
reported  in  [15]  for  the  case  of  R  =  0, 

Both  of.  the  theoretical  velocity  histories  shown  in  Figure  11  are 
in  reasonable  agreement  with  the  primary  characteristics  of  the>  experi¬ 
mental  histories.  The  influence  of  the  added  mass  can  be  seen  by 
comparing  the  upper  and  lower  histories  of  Figure  11.  For  the  case  of 
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added  mass,  the  velocity  generally  is  at  a  higher  level  and  does  not 
oscillate  to  the  degree  observed  for  the  case  of  the  uniform  cylinder 
without,  added  mass.  This  can  be  attributed  to  increased  amplitude  in 
the  rigid  body  translation  mode  when  the  added  mass  is  present  and 
reduced  participation  of  the  xextensional  class  oscillatory  modes,  as 
shown  by  the  mode  amplitude  constants  given  in  Table  3. 

The  theoretical  solutions  shown  in  Figure  11  were  obtained  by 
summation  through  the  eighth' mo.de  Oi,  n  =  8).  Solutions  were  also 
obtained  by  summation  through  the  sixth  mode  and  the  difference  between 
the,  sixth  and  eighth  mode  solutions  were  found  to  be  on  t"  order  of 
one  percent i 

3.  Ring  without  Added  Mass 

The  .theory  will  now  be  compared  with  experimental  data  [IB] 
pbtained  from  strain  gages  on,^the.  inside  surface  of  a  6061  aluminum, 

7  inch  outside  diameter,  1/8  inch  thick,  one  inch  wide  ring  impulsively 
loaded  by  a  magnetically-driven  flyer  plate..  The  pertinent  parameters 
are: 

E  =  10  x  106  (psi)- 

a  =  3.4375  (inches) 

m  =  3.1703  x  10~5  --|-C-  ) 

in 

Z  =  0.00011 

R  =  0 

Substituting  equations  (141)  and  (142),  with  appropriate  spatial 
differentiations,  into- equations  (22)  and  (23),,  $he  stress  resultants 
are  obtained, 
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T  .  -  BA  ( 

T  “  ",ot{- 


l  +  Z  . 

smwcsO'  + 


COS  0  sill  Wes|  i 


TJ(jJ 


csO 


csl 


,,2'^  t  |7.(ir-,l)- 1 1  r’l 

,£bi-.7r(1-n  >  L  J 


(148) 


G  =  sv0UAZ 


•cos  nor' 

71 


where 


and 


(n2  -  Cfs|1)  sin  <dfsnt  (n2-CCsn)sin<ocsnt- 
wfsn  K‘r  “  Gfsn)-  +  n”i  wesn  i(,r  “  Gesr^-  +  nr  1 

(,|2-cfsn)“«n<Jfsiit  +  ,(»2 -  W^in  coesnt 
wfsn^n”  -  Gfsn^”  +  ll”J  wcsn  ^ir  “  Gesn)~  *  n”  1 


(149) 


To,  obtain  non-dimensional  stresses  on  the  inside  surface  of  the 
ring,  the  following  expressions  are  employed: 


tfTh  _  T 
'oc  PWqAc 


where 


«T  is  the 
oG  is  the 
li  is:  the 
c  is  the 
1q  is  the 


°Gh  _  6G 
1qc  pIiwqAc 

stress  due  to  thrust 

stress  due  to  bending 

ring  thickness 

material  sound  velocity, 

maximum  impulse  amplitude,  Ig-phwg 
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Thus,  the  total  stress  in  non-dimensional  form  at  any  position 
on  the  inside  surface  of  the  ring  is  given  by 


oli  o-jJi  ach 
l0c  IqC  IqC 


(150) 


It  is  also  desired  to  use  non-dimensional  time  as  follows: 

ct 

7  ~  a 

Equation  (150),  summed  through  eight  modes,  is  shown  in  graphical 
form  on  Figure  12  for  two  spatial  positions,  0=0  and  0  =  ~ 

Shown  also  in  Figure  12  -are  the  experimental  results  from  [16.] .  There 
is  good  agreement  between  theory  and  experiment,  particularly  at  the 
0=0  position.  As  an  added  observation,  the  initial  portion  of  the 
stress  history  at  the  0=0  position  indicates  the  effect  of 
"destructive  interference"  of  the  modes  which  can  be  construed  as  an 
index  of  solution  convergence.  Note  that  there  is  no  significant  stress 
'until  after  t  =  ~  ,  which,  based  on-,  elementary  wave  propagation 

theory,  is  the  time  required  for  a  disturbance  to  travel  1/4  of  the 
circumference  of  the  ring. 
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0  =  J  ,  INSIDE  RING 

Figure  12  -  Theoretical  and  Experimental  [16]  Stress  Histories 

R  =  0,  Z  =  0.00011 


D.  Discussion 


In  the  above  analytical  solutions  the  loading  given  by  equations 
(137)  is  assumed  to  occur  instantaneously  over  half  the  circumference. 
Such  conditions  are  never  physically  attainable  in  experiments  since 
the  loading  always  occurs  over  a  finite  period  of  time.  However,  if 
the  duration  of  load  application  is  small  in  comparison  with  the  period 
of  the  lowest  natural  frequency  of  the  system,  then  an  instantaneous  or 
impulse  type  load  condition  is  usually  considered  to  be  a  valid 
approximation  to  the  physical  problem.  In  Case  3  above,  where  a  flyer 
plate  impacts  a  ring  in  air,  an  assumed  instantaneous  loading  is  valid 
as  evidenced  by  the  agreement  between  theory  and  experiment  shown  in 
Figure  12.  For  the  Navy  experiments,  Cases  1  ana  2,  an  impulse  type 
loading  seemingly  is  more  of  an  approximation  to  the  actual  conditions 
and  further  explanation  is  necessary.  A  strong,  rapidly,  decaying  shock- 
wave  interacting  with  a  structure  submerged  in  ?  fluid  medium  is  an 
extremely  complex  non-linear  problem  in  the  complete  sens  .  The 
principal  cause  of  non-linearity  is  the  occurrence  of  cavitation  or 
separation  between  the  structure  and  surrounding  fluids  after  the  shock- 
wave  arrives  at  the  structure.  Cavitation  occurs  because  of  the 
inability  of  the  fluid  to  withstand  any  appreciable  negative  pressure 
resulting  from  the  radiated  rarefaction  wave.  It  is.  largely  the 
presence  of  cavitation,  however,  that  permits  the  application  of  an 
impulse  type  solution  to  describe  the  response  of  the  cylinder  for 
early  times.  From  experiments  [21],  it  has  been  observed  that  the 
first  peak  velocity,  ,  occurs  very  early  (on  the  order  of  1/2  msec) 

after  shockwave  arrival  and  can  be  predicted  solely  from  the  shockwave 
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parameters  and  Mieil  ma3s  properties  without  regard  to  shell  stiffness 
characteristics.  During  this  early  time  period  after  shockwave  arrival, 
the  first  disturbance  in  the  shell  has  propagated  along  the  circumference 
in  both  directions  about  the  initial  contact  point  a  distance  nearly 
equal  to  one-quarter  of  the  circumference  of  the  cylinder.  By  the  time 
the  peak  velocity  is  reached  at  the  first  contact  point,  a  spatial 

distribution  of  velocity  is  present  over  about  one-half  the  circum¬ 
ference  of  the  shell,  Tor  analytical  convenience,  it  is  assumed  that 
this  initial  velocity  is  spatially  distributed  in  the  manner  described 
by  equations  (137).  It  is  noted  that  the  time  required  to  develop  this 
spatial  velocity  distribution  is  small-  in  comparison  with  the  periods 
of  the  lowest  natural  frequencies  of  either  the  extensional  or  flexural 
classes  of  vibration.  An  impulse  response  solution  for  this  particular 
problem  involves  another  approximation,  which  is  that  the  influence  of 
the  shockwave  in  the  surrounding  water  is,  not  accounted  for  except  in 
the  determination  of  wQ  .  the  effect  of  the,  shockwave  as  it  engulfs 
the  cylinder  is  considerably  reduced  because  the  faster  moving  dis¬ 
turbance  in  the  shell  creates  a  precursor  or  a  radiated  wave  in  xhe 
water  ahead  of  the  free-water  shockwave.  Thus  the  water  is  disturbed 
ahead  of  the  propagating  shockwave  which  lessens  the  shockwave 
influence  on  the  cylinder.  Two  other  factors  which  tend  to  reduce  the 
effect  of  the  shockwave  are  the  natural  decay  in  amplitude  of  the  wave 
front  as  the  distance  from  the  source  is  increased,  and  the  changing 
angle  of  incidence  of  the  wave  front  with  the  cylinder.  On.  the  far 
side  of  the  cylinder,  or  shadow  zone,  the  wave  is  considerably  diffused 
due  to  diffraction  over  the  cylinder  as  an  obstacle,  and  due  to  the 
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effect  of  the  radiated  wave  as  described  above.  From  experimental  [21.] 
pressure  measurements  taken  near  the  fluid-cylinder  interface,  cavita¬ 
tion  has  been  observed  to  begin  at  the  time  %  is  reached  and 
persist  for  a  duration  of  several  milliseconds.  During  this  early 
period  of  time,  the  cylinder  behaves  as  though  it  were  decoupled  from 
the  surrounding  fluid  medium.  It  is  this  decoupling,  that  permits  the 
early  response  features  of  the  cylinder  to  be  predicted  reasonably  well 
(see  Figure  11)  by  the  theory  described  herein.  In  other  situations 
where  little  or  no  cavitation  occurs  the  present  theory  would  be  less 
applicable,  if  at  all. 
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V.  DISCUSSION 


Technical  aspects  of  the  natural  vibration  and  impulse  response 
solutions  have  been  discussed  as  appropriate  in  previous  chapters.  The 
discussion  in  this  chapter  will  pertain  to  implications  of  the  theoreti¬ 
cal  results  with  respect  to  state  of  the  art,  current  applications  and 
future  research. 

The  solution  of  a  fundamental  problem  has  been  obtained  herein, 
thus  the  state  of  the  art  in  the  area  of  vibrations  of  mass-loaded 
structures  has  been  advanced.  Heretofore,  as  noted  in  the  literature 
survey,  the  only  solutions  available,  other  than  a  simplified  theoreti¬ 
cal  treatment,  v<ere  those  obtained  by  approximate  methods  of  analysis, 
With  the  solution  for  the  natural  vibrations, a  variety  of  initial 
value  and  forced  vibration  problems  can  now  be  approached  in  a  more 
rigorous  manner.  An  application  of  the  theory  to  an  initial  value 
problem  was  demonstrated  in  the  previous  chapter  and  judging  by  the 
recent  increased  interest  shown  in  the  literature,  there  are  several 
other  problem  areas  in  which  the  present  theory  would  be  directly 
applicable. 

The  effect  of  imperfections  on  the  vibrations  of  bodies,  of 
revolution  is  a  problem  area  of  basic  physical  importance.  Preliminary 
experimental  investigations  [1,2,3]  into  the  vibrations  of  imperfect 
bodies  of  revolution  have  shown  that  the  presence  of  a  small  asymmetry 
tended  to  resolve  a  single  natural  frequency  into  two  nearby  values. 

The  larger  the  asymmetry  or  irregularity,  the  farther  apart  and  more 
distinct  the  frequencies.  Resonance  tests  have  revealed  the  existence 
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of  so-called  preferential  planes,  i.e,,  if  the  body  was  excited  in  one 
of  these  planes,  only  one  resonance  peak  occurred,  but  if  excited  at 
another  position  around  the  body,  two  peaks  were  observed  for  the  same 
mode  number.  For  a  body  of  revolution  with  a  mass  imperfection,  the 
results  of  this  investigation  tend  to  support  such  experimental 
phenomena  since  the  symmetrical  and  anti-symmetrical  branches  of  the 
same  mode  were  found  to  have  different  natural  frequencies,  A 
systematic  combined  theoretical  and  experimental  study  of  such  phenomena 
would  be  a  contribution  as  apparently  no  such  study  has  been  reported  in 
the  literature.  The  governing  equations  and  their  solution  given 
herein  could  be  utilized  in  a  theoretical  study  of  the  effects  of 
combined  mass  and  elastic  imperfections  for  the  two  dimensional 
problem. 

With  regard  to  the  problem  of  an  underwater  shock  wave  interacting 
with  a  cylinder,  some  further  effort  is  needed  to  establish  the  range 
of  applicability  of  an  impulse  type  solution  as  presented  herein.  If 
large  non-linear  cavitation  effects  are  not  present,  then  linear 
acoustic  theory,  recently  extended  by-  Huang  [22],  would  be  more 
appropriate  to  this  problem.  To  determine  their  range  of  validity, 
parametric  studies  should  be  made  with  .both  theories,  impulse  response 
and  acoustic  interaction,  for  comparison  with  various  experimental 
results. 

Concerning  current  applications  of  the  theory,  perhaps  the  most 
important  aspect  from  an  engineering  standpoint  is  that  the  solution 
can  serve  as  a  base  for  comparison  with  approximate  analytical  methods, 
such  as  the  finite  element.  Such  methods  offer  a  means  for  obtaining 
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practical  solutions  of  more  complex  problems.  For  example,  an  analyti¬ 
cal  study  of  the  vibrations  of  a  cylinder  containing  several  lumped 
masses  and  attached  mechanical  systems  would  be  intractable  except 
by  an  approximate  approach.  In  this  connection,  it  is  noted  that  in  an 
approximate  energy  approach  to  such  problems,  improved  trial  functions 
would  be  available  since  the  general  character  of  the  solution  would 
be  known  from  the  present  theory. 

As  the  first  extension  of  the  theory  for  future  research,  it 
would  be  desirable  to  include  rotatory  inertia  of  the  added  mass.  In 
principle,  this  additional  feature  would  offer  little  analytical 
difficulty,  and  the  solution  would  be  more  realistic,  particularly  for 
anti-symmetrical  branch  motions.  Another  future  research  effort  that 
should  be  considered  is  extension  of  the  theory*  to  include  rotatory- 
inertia  and  shear  deformation  effects  .in  the  shell.  Extension  of  the 
theory  in  this  manner,  whether  done  approximately  or  by  using  higher 
order  shell  equations,  would  involve  considerable  analytical  effort; 
however,  the  inclusion  of  such  effects  would  offer  a  more  accurate 
account  of  the  coupling  of  the  extensional  and  flexural  motions  as 
shown  in  Figure  7. 

To  obtain  a  better  physical  representation  of  a  large  added  mass, 
the  theory  should  be  extended  to  include  the  added  mass  distributed 
over  a  finite  length  of  the  circumference,  or  a  finite  length  rigid 
inclusion,  rather  than  concentrated  at  a  point  on  the  circumference. 

If  considering  the  problem  of  a  mechanical  system  attached  to  the 
cylinder,  then  a  more  accurate  treatment  would  be  to  model  both  mass 
and  elastic  properties  of  the  system,  instead  of  just  the  mass  alone. 
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For  each  of  the  problems  mentioned  above,  th,s  governing,  equations  and 
their  solution,  equations  (10)  and  (11),  would  be  applicable,  although 
the  eigenvalue  problem  would  change  due  to  differences  in  boundary 
conditions. 

The  theory  could  perhaps  be  extended  to  other  more  complex 
problems,  such  as  a  finite  cylindrical  shell  with  an  attached  mechanical 
system,  to  mention  one,  but  an  extensive  preliminary  effort  would  have 
to  be  made  in  each  case  to  determine  if  a  rigorous  analytical  treatment 

would  be  tractable. 


VI .  SUMMARY 


the  complete  plane  strain  solution  was  obtained  for  the  natural 
vibrations  of  a  thin  circular  cylindrical  shell  containing  an  added 
concentrated  line  mass.  For  a  particular  mode,  four  distinct  solution 
states  were  found  to  exist,  a  symmetrical  and  anti-symmetrical  branch 
for  each  class,  flexural  and  extensional.  Significant  features 
revealed  by  this  investigation  were  the  difference  in  frequency  and 
mode  shape  of  each  solution  state  and  the  presence  of  coupling  between 
the  flexural  and  extensional  classes,  particularly  noticeable  in  the 
extensional  class  mode  shapes.  The  solution  for  the  case  of  a  uniform 
cylinder  was  obtained  from  the  general  solution  by  taking  the  added 
mass  to  be  zero.  For  the  uniform  cylinder,  the  frequencies  of  each 
branch  of  the  same  class  were,  as  expected,  the  same.  The  influence  of 
the  added  mass  was  determined  by  comparing  numerical  solutions  with 
and  without  the  mass.  It  was  found  that  with  the  added  mass  the 
natural  frequency  was  reduced  for  each  solution  state  with  greater 
reduction  occurring  in  the  flexural  class  than  in  the  extensional  classt 
for  a  particular  mode.  The  results  of  an  experimental  natural  vibra¬ 
tion  study  [8]  were  found  to  be  in  general  agreement  with  the  theory. 

Using  the  natural  vibration  solution,  the  general  impulse 
response  solution  for  arbitrary  initial  conditions  was  obtained  by 
normal  mode  theory.  From  the  general  solution,  numerical  solutions 
for  specified  initial  'conditions  were  obtained  and  compared  with  the 
results  of  several  experimental  studies  [15,16,20],  In  all  cases,  the 
theory  was  found  to  agree  reasonably  well  with  experiments  The  major 
influence  of  the  added  mass  on  the  impluse  response  was  found  to  be  an 


increased  participation  of  the  rigid  body  translation  and  flexural 
class  modes,  and  decreased  participation  of  the  extensional  class 
oscillatory  modes. 

In  the  discussion  of  the  previous  chapter,  several  direct  practi¬ 
cal  applications  of  the  theory  were  mentioned.  The  theory  would  be 
applicable  in  basic  studies  of  the  vibrations  of  bodies  of  revolution 
containing  certain  imperfections,  in  predicting  the  early  response  of 
submerged  cylinders  to  shock  waves  in  the  presence  of  cavitation,  and 
as  a  base  for  comparison  with  approximate  analytical  methods.  Dis¬ 
cussed  also  were  areas  of  future  research  for  extension  of  the  theory 
so  that  it  would  be  physically  more  meaningful  in  certain  future 
applications.  In  this  regard,  it  was  pointed  out  that  an  effort 
should  be  made  to  include  rotatory  inertia  of  the  added  mass  and  to 
consider  rotatory  inertia  and  shear  deformation  effects  in  the  shell. 

It  was  also  pointed  out  that  for  some  situations,  more  realistic 
models  of  the  attached  mass  other  than  as  a  concentrated  line  mass 
should  be  considered  as  topics  for  future  research. 
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